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A density dependent relativistic mean-field model is determined to reproduce the components of 
the nucleon self-energy at low densities. This model is used to investigate spinodal instabilities in 
isospin asymmetric nuclear matter at finite temperatures. The inhomogeneous density distributions 
in the spinodal region are investigated through calculations in a cubic Wigner-Seitz cell. Compared 
to results obtained in phenomenological calculations the spinodal region is large, i.e. the spinodal 
region at zero temperature can reach densities above 0.12 fm -3 . The predicted spinodal region is 
s ! ■ concentrated around isospin symmetric nuclear matter and the critical temperature is considerably 

lower than in the previous microscopic based investigation within a non-relativistic Brueckner- 
. Hartree-Fock approach. 
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I. INTRODUCTION 



The properties of isospin asymmetric nuclear matter systems like their instabilities and phase transitions are 
presently a topic of great interest [U-Q ■ At either high energy or high baryon density one explores the phase transition 
to quark matter, i.e. the hadrons dissolve into a phase of de-confined quarks. Below a critical temperature or density, 
the homogeneous fluid consisting of nucleons may condense into nuclei. This transition can be regarded to be of the 
liquid-gas type. Heavy ion collisions are well suited to investigate these transitions, since the multi-fragmentation 
phenomenon occurs in these experiments. In the multi-fragmentation phenomenon, highly excited composed nuclei 
are formed in a gas of evaporated particles. Therefore, this multi-fragmentation phenomenon can be interpreted as 
the coexistence of a liquid and a gas phase. This liquid-gas phase transition is important in the description of the 
—1 ■ nucleation processes. Furthermore, it plays an important role in the collapse of supernovae into neutron stars. 
^ ' The hadrons involved in this transition can be either neutrons or protons. Therefore, one should consider nuclear 
. matter as composed of two different fluids. These phase transitions in nuclear matter are related to the thermodynamic 
instability regions, which are limited by so-called spinodals. This spinodal instability is dominated by mechanical or 
isoscalar density fluctuations which lead to a liquid-gas phase separation with restoration of the isospin symmetry in 
the dense phase. This phenomenon is known as the isospin distillation or fragmentation effect 0], i.e. the formation 
of large droplets of high density symmetric nuclear matter within a neutron gas with only a very small fraction of 
protons. 

In the past systematic analysis of the stability conditions of asymmetric nuclear matter against the liquid-gas phase 
transition have been performed using different methods. These approaches can roughly be divided into phenomeno- 
logical methods and microscopic based methods. Phenomenological studies of spinodal instabilities involve mean-field 
calculations with effective forces of the Skyrme or Gogny type JlL |9|Jll| and relativistic mean- field calculations using 
constant and density dependent coupling parameters [If @[, 0, Il3j|. Some of these models predict that the total 
density at which spinodal instability appears will decrease if the asymmetry increases, whereas for other models it 
will increase up to large asymmetry and will finally decrease. However, a common feature of these models is that 
the spinodal region is in principle situated at densities below ps < 0.1 fm~ 3 . The microscopic based studies are a 
calculation with low-density functionals || based on the Dirac-Brueckner-Hartree-Fock (DBHF) approach [l4| with 



a realistic Bonn A potential [15[ and a calculation within a microscopic Brueckner-Hartree-Fock (BHF) approach ^ 
using the realistic Argonne V18 potential [l6| plus a three-body force of Urbana type. Both calculations predict larger 
spinodal regions compared to those obtained in the phenomenological calculations, i.e. the spinodal region can reach 
densities above 0.12 fm -3 close to symmetric nuclear matter. The critical temperature in the BHF calculation is 
T c = 17.5 MeV, whereas finite temperature effects are not considered in the work of Ref. 

In the present work, we construct a density dependent relativistic mean-field (DDRMF) model, i.e. a density 
dependent relativistic Hartree (DDRH) model, which is based on microscopic DBHF calculations [l7| using the realistic 
Bonn A potential [l5| • Apart from the application of the results of this DBHF calculation to spinodal instabilities [|| , 
it has also already been used for other applications such as investigations of the non-homogeneous structures in the 
neutron star crust [l8[, investigations of hybrid stars [Hj], and studies of nucleon-nucleus scattering [2(J. In the 
previous DBHF based investigation of the spinodal instability Q , a functional is constructed on a density expansion 
of the relativistic mean-field theory with fixed coupling constants. However, in the present work the approximation 
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of an expansion is not made and the effective coupling functions, which are obtained from the DBHF self-energy 
components, are density dependent. Special attention is devoted to reproduce the DBHF self-energy also at very low 
densities, i.e. ps < 0.10/m -3 . In the framework of this DDRH theory based on a DBHF approach, the equation of 
state at zero and finite temperature can then be determined from an iteration procedure until convergence is reached. 
Within this model, the temperature effects and the nature of spinodal instabilities are investigated. In addition to 
these investigations in infinite nuclear matter calculations, the density distribution in the spinodal region is studied 
in cubic Wigner-Seitz (WS) cell calculations. 

The paper is organized as follows. The density dependent relativistic Hartree (DDRH) theory is treated in Sec. [TT1 
The parameterization for the DDRH theory based on the DBHF results are discussed in Sec. Mil Section lTVl is devoted 
to the theory of the spinodal instabilities. The results for the spinodal instabilities are presented and discussed in 
Sec. [V] Finally, we end with a summary and the conclusion in Sec. IVII 



II. DENSITY DEPENDENT RELATIVISTIC HARTREE THEORY. 

A Lagrangian density of an interacting many-particle system consisting of nucleons and mesons is the starting 
point of a DDRH theory. The Lagrangian density of the DDRH theory used for our investigation includes as well 
the isoscalar mesons a and u) as the isovector mesons S and p. Therefore, the Lagrangian density has the following 
parts: the free baryon Lagrangian density Cb, the free meson Lagrangian density Cm-, and the interaction Lagrangian 
density £ int : 



which can be written as 



£b + £m + Ant, (1) 



£ B = *(n^-M)*, (2) 



£m = I ( W"** - m?*?) - | £ (\F {K) , V Fft KA^A^) , (3) 



Ant = -g<j^<5>„^ - 35*t* 5 * - <^*7ju^(L)* - 9p^r^A^. (4) 

In the above Lagrangian density the field strength tensor for the vector mesons is defined as .F( re w = d IJi Ai K \ v —d v Ai K \ li . 
The nucleon field is ^> and the nucleon rest mass is M. The meson fields are denoted by A^), and A^ p y 

Moreover, the bold symbols denote vectors in the isospin space. The mesons couple to the nucleons with the strength 
of the coupling constants g a , gs, g^, and g p . The rest masses of these mesons are m CT , m^, ms, and ra p . 

By minimizing the action for variations of the fields ^ included in the Lagrangian density of Eq. ([1} , one obtains 
the field equations for the nucleons, 

5 dt J d 3 x£($(x),d^(x),t) =0. (5) 

Finally the following field equation can be derived 

d dC dC dC 5p 



dx» d(d^) 9* dp 



0. (6) 



Without density dependent meson-baryon vertices the third term vanishes in Eq. ([6]) and the normal Euler-Lagrange 
field equation is obtained. From the evaluation of Eq. © , we obtain the Dirac equation, 

(i 7M ^ - M - (S + E< r >7b)) * = 0, (7) 

with E the nucleon self-energy and E^ r ' the so-called rearrangement contribution to the nucleon self-energy. 
The self-energy E for the Lagrangian density in Eqs. ©-(HI) can be written as 



E = g a <$> a + g S T<f>s + 9ul^ u) + g p T^ p A^ p) . 



(8) 
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behavior under Lorentz transformations, 

E(fc) = S s (fc)- 70 S o (fc)+ 7 .kI] v (fc). (9) 
The S s , S OJ and S v components are Lorentz scalar functions. Hence, one can define the following effective quantities 

k* = k(l +&eE v (k)), (fO) 

m *(k) = M + 5ReS s (fc), (ff) 

and 

E*(k) =E(k) + $teX (k). (12) 
In the on-shell case, the effective energy can also be given by 

E* (k) 2 = k* 2 + m* (fc) 2 . (13) 
The contributions to the self-energy of the Lagrangian density presented in Eqs. (|TJ)-(|3J) can be written as 

EmO) = -[ — ) (Ps,n + Ps,p) - [ — ) Y" (ps,i ~ Ps,j) 

and 

s 0ii (fc) = - ( Pn + Pp ) - (js.) x: 

i r . ,.,,2 m * 



and 

Pi = ^ / Mk)k 2 dk. (15) 
The thermal occupation factors in Eqs. (fT4|) and ([TBI are given by 

^ (fc)= l + ex P [(4-^)/T]' i=P '» (16) 

in the "no-sea" approximation. 

The term in Eq. (|7J| is generated by the third term in Eq. ©. It will only give a nonzero contribution, if 
meson-baryon vertices are density dependent. This rearrangement contribution £( r ) can be written as 



+ t i% ^)* + ^ fT7 ^*)- (17) 

Rearrangement contributions are essential to provide a symmetry conserving approach pll [22l |. i.e. an approach 
satisfying energy-momentum conservation and thermodynamic consistency like the Hugenholtz - van Hove theorem. 

III. PARAMETERIZATION 



We have constructed a density dependent relativistic mean-field (DDRMF) model, i.e. a density dependent relativis- 
tic Hartree (DDRH) model. The effective density dependent coupling functions of this DDRH model are determined 
by the requirement that the results of a DBHF calculation at zero temperature using Bonn A [17j for the components 



4 



meson i 


m [MeV] 


at 


bi 


a 


di 


e. 


a 
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14.0089064 


9.91395378 


-2.1599853 


0.280160069 


-13.3718958 




782.6 


13.5288515 


7.86182976 


-1.8871752 


0.244318038 


-9.42372704 


P 


769 


4.90660906 


-14.9182243 


3.93277073 


-0.498567462 


13.1945915 


s 


983 


-8.32733536 


-43.6584053 


10.4067049 


-1.26465869 


50.3025665 



TABLE I: Parameter set for the DDRH model from the DBHF approach in Ref. [17 
400 r 




0.2 0.3 
density [fm 3 ] 



0.5 



FIG. 1: (Color online) Density dependence of the coupling functions of the cr-meson (solid), the cj-meson (dotted), the p-meson 
(dashed), and the 5-meson (dashed-dotted). 



of the nucleon self-energy at the Fermi momentum are reproduced. The rearrangement terms of DDRH are not 
included in the fit, since the DBHF approach has no rearrangement contributions. 
The density dependence of the coupling functions are parameterized by 

9i{pb) = a% + M + CiX 2 + diX 3 + e l yfx, (18) 

with x = pb/po and p = 0.16 fur 3 . The DBHF data at densities of p B = 0.017, 0.100, 0.197, 0.313, and 0.467 fm" 3 
at a proton fraction of Y p = 0.4 are used to obtain these parameters of our DDRH model. The meson masses are 
chosen to be identical to those of the Bonn A potential. All parameters are given in table fl] The density dependence 
of these coupling functions is displayed in Fig. Q] The density dependence is rather weak except at very low densities. 
The strength for the isovector mesons (p and S) is significantly smaller than for the isoscalar mesons (a and to). All 
coupling functions tend to decrease with increasing density. This may reflect the quenching of the DBHF G-matrix 
underlying this fit due to Pauli and dispersive effects in the Bethe-Goldstone equation. 

In Fig. [2] the binding energy of the DDRH model is shown for neutron matter and symmetric nuclear matter. The 
results are in reasonable agreement with the DBHF results. Note, that the DDRH parameterization has solely been 
based on the DBHF self-energies at the corresponding Fermi momenta for a proton fraction Y p — 0.4. No attempt 
has been made to account for the momentum and further asymmetry dependence of the DBHF self-energies, which 
enter the calculated energies for symmetric nuclear matter and neutron matter. Both approaches show a nonlinear 
convergence to zero when the density decreases, which is in qualitative agreement with a study of low-density nuclear 
matter [23j based on a virial expansion that includes protons, neutrons, and a-particle degrees of freedom. In particular 
at densities below p < 0.1fm~ 3 , the DDRH model presented here is in better agreement with the DBHF result than 



those models constructed in Ref. 24 



From table [TT1 one can see that the saturation density p sat is shifted to lower densities and the binding energy 
E/A of the saturation point is weaker in the DDRH model compared to the original DBHF results. The compression 
modulus K is somewhat smaller and symmetry energy a s is in a fairly good agreement with the DBHF results. 
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FIG. 2: (Color online) Binding energy is plotted as a function of the density for the DBHF approach and the DDRH model. 
Left: isospin symmetric nuclear matter. Right: pure neutron matter. 





DDRH 


DBHF 


psat [fm 3 ] 


0.1759 


0.181 


E/A [MeV] 


-14.39 


-16.15 


K [MeV] 


193 


230 


a s [MeV] 


35.03 


34.16 



TABLE II: Saturation properties of nuclear matter for the DDRH model compared to the DBHF results of Ref. [3, fl7j . 
The quantities listed include the saturation density p sa t, the binding energy E/A, the compressibility modulus K, and the 
symmetry energy a s . The symmetry energy a„ has been obtained assuming a quadratic dependence of the EOS on the 
asymmetry parameter /3 — (p n — p P )/(p n + p P )- 



IV. SPINODAL INSTABILITIES 



The stability conditions for isospin asymmetric nuclear matter at finite temperatures are obtained from the free 
energy density T . This free energy density is defined as 

T = £ - TS, (19) 

where £ is the energy density and S denotes the entropy density. The energy density £ is obtained from the energy- 
momentum tensor 

£=(T m ) = -L J fi(k)k 2 dk y/k*+m? 

i—p : n 

where /,•(&) is the thermal occupation probability. The rearrangement term does not appear in Eq. (|20[) . since 
rearrangement does not contribute to the energy density. However, note that the rearrangement does modify other 
quantities such as the pressure. Furthermore, the entropy density S is given by 

5 = -^E / /*( fc ) fc2 dk t^( fc ) Hh(k)) + (1 - /*(*)) In(l - fi{k))\ . (21) 

i— p,n 

The necessary ingredients to analyze the stability of isospin asymmetric nuclear matter against separation into two 
phases are now introduced. The nuclear system will be stable if the free energy of a single phase is lower than the 
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free energy in any two-phase configurations. This stability condition implies that the free energy density is a convex 
function of the neutron and proton densities p n and p p . This means that the curvature matrix, 



d 2 T 



dpidpj 



(22) 



is positive-definite. The requirement that the local curvature is positive demands that both the trace and the deter- 
minant are positive, 



Tr[J" y ] > 0, and Det[J"y] > 0. 



(23) 



The [Fij] is a 2 * 2 symmetric matrix, since one considers a two-fluids system. Therefore, the stability matrix [Tij] 
has two real eigenvalues : 



associated to eigenvectors Sp given by (i ^ j) 



X ± = I ( TvL/V ± 



Tr[^-] 2 -4Det [T i3 \ 



Sp 



A ± — T, 



33 



A — J~i, 



Therefore, the requirements in Eq. (|23j) can equivalently be written as 

Tr[J"y] = A+ + A" > and DetfJy] = A+A~ > 0. 



(24) 



(25) 



(26) 



This means that both eigenvalues should be positive to guarantee the stability of the system. Due to the essentially 
quadratic dependence of the symmetry energy on the asymmetry parameter and a positive increasing symmetry 
energy with the total density, it turns out that only one eigenvalue can become negative [![. In principle A + will 
always be positive in isospin asymmetric nuclear. Only A~ can become negative and the eigenvector associated with 
this negative eigenvalue indicates the direction of the instability. Therefore, if Eq. ([26]) is violated, or equivalently A~ 
is negative for isospin asymmetric nuclear matter, the system will be in the unstable region of a phase transition. 



V. RESULTS AND DISCUSSION 



In the present work, we have investigated the spinodal region in the case of isospin asymmetric nuclear matter at 
finite temperature. The projection of the spinodal contour in the neutron-proton density plane for several temperatures 
ranging from T=0 to T=10 MeV is depicted in Fig. [3l Keep in mind that the isospin symmetry of the nucleon-nucleon 
interaction guarantees the symmetry of the spinodal contours with respect to the interchange of the neutron density 
p n and the proton density p p . 

Let us first consider the results for zero temperature (T= MeV). One finds that the spinodal region extends up to a 
maximum density of 0.127 fm~ 3 . This is considerably larger than the maximum density obtained in phenomenological 
models [lHj, |7J, 191413] . which indicate that the spinodal region is situated at densities below ps < 0.1 fm~ 3 . The 
large spinodal region, however, is in agreement with microscopic calculations 0, @. This difference between the 
phenomenological and the microscopic based models is of course related to the larger saturation densities, which are 
obtained in the microscopic calculations. Note, however, that the present study yields a larger maximal density than 
the non-relativistic BHF calculations of Vidana and Polls [6j , although the saturation density for the present approach 
(po = 0.176 fm~ 3 ) is smaller than the one obtained in the BHF calculation (po — 0.182 fm~ 3 ). 

The shape of the spinodal region for T=0 projected in the neutron-proton density plane, is similar to the shapes 
obtained in calculations before. In detail, however, there is a difference: The region of instability is more confined to 
isospin symmetric matter. 

Also the dependence of the region of instability on temperature in Fig. [3] is very similar to the one observed 
in previous calculations flUH 13. 171 J9Tfl3j . The spinodal region shrinks with increasing temperature, until the critical 
temperature T c is reached. In the present approach no spinodal region can be obtained above T — 11 MeV. Therefore, 
the critical temperature is considerably lower than the value T c — 17.5 MeV obtained in the other microscopic based 
calculation Q or the values for T c reported in the phenomenological studies It is worth mentioning that the 
spinodal region obtained here at temperatures close to T c occurs at slightly higher densities than in the other studies. 

It has been shown that the spinodal instabilities in isospin asymmetric nuclear appear as a mixture of isoscalar and 
isovector instabilities [H, Hj]. Isoscalar or mechanical instabilities are related to density fluctuations. Pure isoscalar 
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FIG. 3: (Color online) The projection of the spinodal contour in the neutron-proton density plane at several temperatures. 




FIG. 4: (Color online) Ratio Spp /5pn as a function of asymmetry parameter ft — (p„ — - p p ) / (p„ + p p ) for several temperatures. 
The left panel shows the results for a fixed density of ps = 0.03 fm -3 , whereas in the right panel the results for a fixed density 
of pB = 0.06 fm -3 are presented. The ratio p v j 'p n is given by the black dotted line. 



instabilities conserve the proton fraction and fulfill the relation 5p~ /5p~ = 1. On the other hand, pure isovector 
or chemical instabilities are related to fluctuations in the proton fraction at fixed density and fulfil the relation 
Sp~ /5p~ = — 1. Therefore, the instabilities will be predominantly of the isoscalar nature, if the neutrons and protons 
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FIG. 5: (Color online) Neutron density distributions at T = MeV resulting from box calculations for isospin symmetric 
nuclear matter as a function of the Cartesian coordinates x,y, and z. The panels refer to densities of p — 0.005 fm~ 3 , 0.05 
fm -3 , 0.10 far 3 , and 0.125 fm -3 . 



move in phase 5p p /6p n > 0, or of the isovector nature, if they move out of phase Sp p /5p n < 0. In Fig. U the 
results obtained from our model for this ratio 5p~ /6p~ is plotted as a function of the asymmetry parameter j3 — 
(fin ~ Pp)/{Pn + Pp) for several temperatures ranging from T — up to T — 8 MeV and for two different densities. It 
can be observed from Fig. 2] that the ratio 8p~ /Sp~ is positive in all considered cases. This result implies that the 
spinodal instability is dominated by mechanical or isoscalar density fluctuations, which is in agreement with previous 
investigations [H, 0, 0, HJ. Since Sp~ /5p~ is larger than the ratio p p /p n presented by the black dotted line in Fig. |4j 
the dense liquid phase of the system is driven towards a more isospin symmetric composition and the gas phase will 
be more isospin asymmetric. These effects are less strong for higher temperatures. The liquid phase is in general also 
less strongly driven towards a isospin symmetric composition, when the density or the asymmetry parameter j3 
increase. 

Up to this point we have determined the spinodal zone from the evaluation of the curvature matrix [J^ij] in (|22p 
for the free energy of the homogeneous system of nuclear matter. Having identified the regions of instability we will 
now try to verify these instabilities and determine the density profile in this area. 

For that purpose we consider the very same DDRH and perform variational calculations in a cubic Wigner-Seitz 
(WS) cell. The size and shape of the WS cell should be determined in a variational manner, however, in a kind 
of explorative study we have chosen the length of the WS cell to be 20 fm in each direction with the origin of the 
coordinate system in the center of the box. Also we have restricted the study to isospin symmetric systems. In order 
to be consistent with our infinite nuclear matter calculations, the Coulomb energy is not included in the WS cell 
calculations. Therefore the proton density distributions are identical to the corresponding neutron ones displayed in 
Fig. El 

The transition from isolated nuclei to a phase of homogeneous baryon matter can be observed in the various panels 
Fig. [5]for zero temperature. The left top panel shows the density distribution at a density of p = 0.005 fm -3 , i.e. only 
20 protons and 20 neutrons ( 40 Ca) are present in the cubic WS cell. At these low densities, isolated nuclei will be 
present. The density profiles are identical in all three Cartesian directions, which corresponds to crystal of spherical 
nuclei located at the centers of the WS cell. 

However, at higher densities this is not the case anymore. At a density of p = 0.05 fm~ 3 , the density profile in 
z-direction deviates a little from that in the x and y direction. Furthermore, one has a high-density region at the 
center of the cubic WS cell and a low-density region along the connection lines to the nearest neighbor at the boundary 
of the box. 

However, the density distribution is not as simple as suggested by Fig. [5] This can be seen from the complete 
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density distribution in the z = plane as displayed in Fig. [6] Apart from the low density regions at the boundary of 
the box along the connection lines between the centers of the WS cells, one has also regions in which the density is 
almost zero. Therefore, the geometrical structure at this density can be described by large quasi-nuclei located in the 
centers of the WS cell, which are connected to each other by low-density bridges. 

At even higher densities, i.e. around p = 0.10 fm~ 3 , the opposite density distribution can be observed, i.e isolated 
low-density regions completely surrounded by high-density matter. At a density of p = 0.125 fm -3 , so just a little bit 
below the critical density (p c = 0.127 fm -3 ) predicted from the study of the curvature matrix in homogeneous one 
obtains a grid structure with small fluctuations, which resembles homogeneous nuclear matter. 

Therefore, within the numerical accuracy we regard this result to be in agreement with the result of the infinite 
nuclear matter calculation in Fig. [31 

Performing DDRH calculations in the same kind of WS cell at the finite temperature T = 8 MeV, we obtain 
the neutron density distributions displayed in Fig. [7J It is clear that finite temperature effects tend to dissolve the 
geometrical structures observed for T = MeV. The structures at T = 8 MeV are less pronounced than those for T 
= MeV in Fig. [6] and the transition density is lower at T = 8 MeV. These results are in agreement with the results 
presented in Fig. [3jof the infinite nuclear matter calculations. 



VI. SUMMARY AND CONCLUSION 



In the present work, we construct a Density Dependent Relativistic Mean Field model (DDRH) to parameterize the 
results of microscopic Dirac Brueckner Hartree Fock (DBHF) calculations [TtJ using the realistic Bonn A potential [lj| . 
Special attention is devoted to reproduce the DBHF self-energy in particular at densities below the saturation density. 
The equation of state at zero and finite temperature can then be determined in the framework of this DDRH theory for 
homogeneous matter as well as inhomogeneous structures. The spinodal instabilities of homogeneous matter including 
their temperature effects are investigated within this model. 

Our investigations confirm the main results of earlier studies based on phenomenological 0, Q as well as realistic 
models of the NN interaction [fj, that the spinodal instability is dominated by mechanical or isoscalar density fluc- 
tuations which lead to a liquid-gas phase separation with restoration of the isospin symmetry in the liquid phase. 




FIG. 6: (Color online) Neutron density distribution at T = MeV in the z = plane resulting from box calculations for isospin 
symmetric nuclear matter. The panels refer to densities of p = 0.005 fm -3 , 0.05 fm -3 , 0.10 fm -3 , and 0.125 fm -3 . 
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p= Q.Q5 fra 




FIG. 7: Neutron density distribution at T = 8 MeV in the z = plane resulting from box calculations for isospin symmetric 
nuclear matter. The panels refer to densities of p = 0.05 fm -3 , 0.10 fan -3 , and 0.125 fan -3 . 



However, this restoration is less strong for higher temperatures. 

Furthermore, the DDRH model predicts some details, which are different from earlier studies: A large critical 
density at zero temperature, a smaller critical temperature T c and details in the shape of the spinodal region. At 
zero temperature the spinodal region reaches densities above 0.12 fm -3 . This is different from the predictions of 
phenomenological calculations [lHj, LD, l9Ml3| , but in line with the microscopic studies of Vidana and Polls [6] . It 
can to some extent be attributed to the differences in the predicted saturation point of nuclear matter. The present 
DDRH approach yields a smaller critical temperature (T c less than 11 MeV) and higher densities for the spinodal 
region at temperatures just below T c . The predicted spinodal region obtained from our calculation is more strongly 
concentrated around isospin symmetric nuclear matter. 

The studies of the spinodal regions which are based on the curvature of the free energy of homogeneous matter 
are supplemented by variational calculations allowing for density fluctuations explicitly. Calculations employing the 
same DDRH model are performed in a cubic Wigner Seitz cell by minimizing the total energy to obtain density 
distributions in isospin symmetric nuclear matter. At a density of 0.125 fm -3 , a geometrical structure resembling 
homogeneous nuclear matter can be found. At lower densities, non-homogeneous nuclear matter structures occur. 
Finite temperature effects tend to dissolve the geometrical structures observed, i.e. the structures at finite temperature 
are less pronounced and the transition density to homogeneous nuclear matter is lower. 
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